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(N ■ ABSTRACT 

■ We present an analysis of the X-ray spectrum of the Galactic shell-type SNR 
(3 ! G347.3-0.5 (RX J1713. 7-3946). This SNR is a member of a growing class of SNRs 

. which are dynamically young, shell-type sources that emit non-thermal X-rays from 
specific regions on their outer shells. By performing a joint spectral analysis of data 

> ; from observations made of G347.3-0.5 using the ROSAT PSPC, the ASCA GIS and 

\^ I the RXTE PGA, we have fit the spectra of particular regions of this SNR (including the 

' bright northwestern and southwestern rims, the northeast rim and the interior diffuse 

■ emission) over the approximate energy range of 0.5 through 30 keV. We find that fits to 
. the spectra of this SNR over this energy range using the SRCUT model were superior to 



o 



a simple power law model or the SRESC model. We find that the inclusion of a thermal 
model with the SRCUT model helps to improve the fit to the observed X-ray spectrum: 
^ , this represents the first detection of thermal X-ray emission from G347.3— 0.5. Thermal 

|H \ emission appears to be more clearly associated with the diffuse emission in the interior 

^ ' of the SNR than with the bright X-ray emitting rims. A weak emission feature seen 

near 6.4 keV in the RXTE PGA spectrum most likely originates from diffuse X-ray 
emission from the surrounding Galactic Ridge rather than from G347.3— 0.5 itself. We 
have analyzed our RXTE PGA data to search for pulsations from a recently discov- 
ered radio pulsar (PSR J1713— 3949) which may be associated with G347.3— 0.5, and 
we do not detect any X-ray pulsations at the measured radio period of 392 ms. Using 
the best-fit parameters obtained from the SRCUT model, we estimate the maximum 
energy of cosmic-ray electrons accelerated by the rims of G347.3— 0.5 to be 19-25 TeV 
(assuming a magnetic field strength of = lO^G), consistent with the results of Ellison 
et al. (2001). We present a broadband (radio to 7-ray) photon energy-flux spectrum for 
the northwestern rim of G347.3— 0.5, where we have fit the spectrum using a more so- 
phisticated synchrotron- inverse Gompton model with a variable magnetic field strength. 
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Our fit derived from tfiis model yields a maximum energy of only S.Stii TeV for the 
accelerated cosmic-ray electrons and a much greater magnetic field strength of ISO^gg'' 
fiG: however, our derived ratio of volumes for TcV emission and X-ray emission based 
on this fit - Vtcv/^ X-ray ~ 1000 - is too large to be physically acceptable. We argue 
that neither non-thermal bremsstrahlung nor neutral pion paritcle decay can adequately 
explain the TeV emission from this rim, and therefore the physical process responsible 
for this emission at this site is currently uncertain. Finally, we compare the gross prop- 
erties of G347.3— 0.5 with other SNRs known to possess X-ray spectra dominated by 
non-thermal emission. 

Subject headings: acceleration of particles, cosmic rays, supernova remnants, X-rays: in- 
dividual (G347.3-0.5 (RX J1713.7-3946), Galactic Ridge, IWGA J1713.4-3949, PSR 
J1713-3949) 



1. Introduction 

Supernova remnants (SNRs) are considered the leading candidates for the accelerators of 
cosmic-rays to at least the "knee" of the cosmic-ray energy spectrum (E « 3000 TeV). While 
the relationship between SNRs and cosmic-ray acceleration has received extensive attention in the 
literature (Bell 1978a,b; Axford 1981; Blandford k Eichler 1987; Jones & Elhson 1991; Malkov & 
Drury 2001), most of this discussion has been theoretical, with only a limited amount of observa- 
tional data to consider. One major unresolved issue related to cosmic-ray acceleration by SNRs is 
the maximum energy E cutoff attained by cosmic-ray particles accelerated by the shock front: for 
samples of young shell-type X-ray emitting SNRs located in the Galaxy and the Large Magellanic 
Cloud, Reynolds & Keohane (1999) and Hcndrick & Reynolds (2001) estimated values for E cutoff 
of < 200 TeV and < 80 TeV, respectively, both well below the cosmic-ray energy spectrum "knee." 
Progress on this issue has been limited due to the lack of both X-ray and TeV observations of SNRs 
with the angular and spectral resolution needed to probe the acceleration process. 

Recently, however, high spatial resolution observations of SNRs at X-ray energies have become 
feasible for the first time using such X-ray observatories as ROSAT, ASCA, Chandra and XMM. 
These high resolution observations yield spatially-resolved spectroscopy of X-ray emitting features 
within the SNR, particularly the X-ray luminous rims that are associated with the expanding shock 
fronts of SNRs. By performing a detailed spectral analysis of X-ray emission from the rims of SNRs, 
corresponding to shock fronts where interstellar or circumstellar material is being swept up, the 
cosmic-ray acceleration process may be probed. Initial insights on cosmic-ray acceleration by SNRs 
have been gleaned from ROSAT, ASCA and RXTE observations of the high-energy X-ray emission 
from such SNRs as Gas A and SN 1006 (Willingalc et al. 1996; Allen et al. 1997; Keohane 1998; 
Dyer et al. 2001; Allen, Petre & Gotthelf 2001), but a thorough multi-wavelength study of SNRs 
is required in order to obtain a detailed understanding of the cosmic-ray acceleration process. 
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In this paper, we present our analysis of the non-thermal X-ray emission from the Galactic 
SNR G347.3-0.5 (also known as RX J1713. 7-3949). First discovered during the Rontgensatellit 
(ROSAT) All Sky Survey (Pfeffcrmann & Aschcnbach 1996), this SNR has attracted a considerable 
amount of interest from the astronomical community. G347.3— 0.5 is one of the brightest Galactic 
X-ray SNRs known with an X-ray flux density of 4.4 x 10~^° erg cm~^ s~-^ (Pfeffermann Sz Aschen- 
bach 1996). Subsequent studies of this SNR based on ROSAT and ASCA observations (Koyama et 
al. 1997; Slane et al. 1999) revealed that most of its X-ray emission is non-thermal and most likely 
produced by the synchrotron process. Unlike the X-ray spectra typically seen from most other 
SNRs, no thermal line emission was detected from any portion of G347.3— 0.5 by these previous 
studies. The detection of non-thermal X-ray emission from this SNR has made G347.3— 0.5 one 
of the prime examples of an SNR acting as an accelerator of cosmic-ray particles. In addition, 
Muraishi et al. (2000) detected emission from TeV photons located on the northwestern rim of 
G347.3— 0.5: those authors argued that these photons were produced by the invcrse-Compton scat- 
tering process, with cosmic microwave background photons scattering off high-energy cosmic-ray 
electrons accelerated along the shock front of this SNR. This result makes G347.3— 0.5 only the 
third shell-type SNR besides SN 1006 (Koyama et al. 1995) and Gas A (Aharonian et al. 2001) 
reported to produce TeV emission. A detailed analysis of the evolutionary state of G347.3— 0.5 
using radio. X-ray and TeV observations of the bright northwestern rim was given by Ellison et 
al. (2001). Those authors also estimated the maximum energy of cosmic-ray protons accelerated 
by G347.3— 0.5 to be approximately 20 TeV. Butt et al. (2001) argued that this SNR is acceler- 
ating cosmic-ray nuclei as well as cosmic-ray electrons based on positional coincidence between 
the northeastern rim of G347.3-0.5 and the 7-ray source 3EG J1714-3857 (Hartman et al. 1999), 
which is positionally coincident with a molecular cloud: this cloud is assumed to be interacting with 
G347.3— 0.5. Gosmic ray acceleration by this SNR was also discussed by Enomoto et al. (2002), who 
presented a TeV spectrum for the northwestern rim of G347.3— 0.5 and argued that this spectrum is 
best fit by the decay of neutral pion (vr*^) particles rather than an inverse Gompton spectrum of up- 
scattered cosmic background photons. These pions are presumed to be produced by the collision of 
cosmic-ray protons (assumed to be accelerated by G347.3— 0.5) with protons in the molecular cloud, 
as described by Enomoto et al. (2002). This result is in part based on the assumption that 3EG 
J1714— 3857 is associated with G347.4-0.5 and, if correct, would represent significant observational 
evidence for cosmic-ray proton acceleration by an SNR (Aharonian 2002). However, Reimer & Pohl 
(2002) and Butt et al. (2002) have disputed the association with the 7-ray source and instead claim 
that, while G347.3— 0.5 is very likely an accelerator of cosmic-ray electrons, no convincing evidence 
yet exists that it accelerates hadrons as well. Other recent work on G347.3— 0.5 includes model- 
ing of the spectrum of AX J1714.1-3912, a hard X-ray source associated with 3EG J1714-3857 
(Uchiyama et al. 2002), the discovery of fine-structure filaments associated with the northwestern 
rim of G347.3— 0.5 (Uchiyama et al. 2003) and finally, a discussion of the flux of high energy neu- 
trinos from the northwestern rim of this SNR, assuming that the TeV spectrum is produced by the 
decay of neutral pions, which yield neutrinos in the decay process (Alvarez-Muiiiz & Halzen 2002). 
In this paper, we present a spatially-resolved broadband spectral analysis of G347.3— 0.5 in order to 
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analyze the high energy emission from this SNR and probe cosmic-ray acceleration by this source. 
Like many Galactic SNRs, the distance to G347.3— 0.5 and its age are not well known: early studies 
of this SNR suggested a low distance of about 1 kpc assuming G347.3— 0.5 is in the Sedov stage of 
evolution (Pfeffermann &: Aschenbach 1996) or that it is associated with the historical supernova of 
393 A.D. (Wang, Qu h Chen 1997). In addition, Koyama et al. (1997) derived a distance of 1 kpc 
by comparing their measurement of the column density toward this source with the total column 
density in the direction of the Galactic Center. Such a low distance suggests that G347.3— 0.5 is a 
young SNR with an age of about 2000 years. However, a subsequent analysis by Slane et al. (1999) 
instead favored a larger distance of 6 kpc based on a postulated association between G347.3— 0.5, 
an HII region and a complex of three molecular clouds: such a distance would imply a much greater 
age of about 10^ years. We will adopt a distance to G347.3— 0.5 of 6 kpc for this paper in order to 
be consistent with the analyses presented by Slane et al. (1999) and Ellison et al. (2001) as well as 
the putative association between this SNR and the recently discovered pulsar J1713-3949, believed 
to lie at a distance of 5.0±0.2 kpc based on measurements of its dispersion measure (Crawford et 
al. 2002). 

The organization of this paper is as follows: in Section 2, we describe the observations of 
G347.3-0.5 using ROSAT (Section 2.1), ASCA (Section 2.2) and RXTE (Section 2.3). While 
the ROSAT and ASCA observations have been presented previously, the RXTE observations are 
presented here for the first time. The spectral fitting and analysis for the datasets from all of 
the observations are described in Section 3, including descriptions of the models used to fit the 
non-thermal and thermal components to the X-ray emission. The results of this spectral fitting are 
presented in Section 4: we discuss in turn the maximum energies of electrons shock-accelerated by 
G347.3— 0.5 (Section 4.1), our search for thermal X-ray emission from different regions of this SNR 
(Section 4.2), and the nature of the weak emission line found in the RXTE PCA spectrum of this 
source (Section 4.3). In our RXTE PCA data, we have searched for pulsed X-ray emission from 
the recently detected radio pulsar which is possibly associated with the discrete X-ray source seen 
at the center of G347.3— 0.5: we describe this search in Section 4.4. A study of the broad-band 
spectrum of the northwestern rim of G347.3— 0.5 from radio to 7-ray wavelengths is presented in 
Section 5. Comparisons between G347.3— 0.5 and two similar SNRs (SN 1006 and G266.2— 1.2) are 
presented in Section 6, and lastly the conclusions are given in Section 7. 



2. Observations and Data Reduction 

In this Section we discuss the X-ray observations that were used in this analysis, including 
the instruments used and the data reduction processes. These observations include archived data 
from observations of G347.3— 0.5 that were made using the Position Sensitive Proprtion Counter 
(PSPCB, hereafter referred to as just PSPC) of ROSAT and the Gas Imaging Spectrometer (CIS) 
of ASCA. We include our observations of this SNR using the Proportional Counting Array (PCA) of 
the Rossi X-ray Timing Explorer [RXTE). These data sample the X-ray emission from G347.3— 0.5 
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over nearly two full decades of energy (0.5 keV to 30 keV). Details of the observations for each 
instrument are provided in Table 1. 

2.1. ROSAT PSPC Observations and Data Reduction 

ROSAT observed G347.3— 0.5 at a location near the center of the remnant using the PSPC, a 
multiple-wire proportional counter sensitive to photons that have energies between approximately 
0.1 and 2.0 keV. The energy resolution AE/E is 0.43 at 0.93 keV and the field of view of the 
telescope is 2°. For our off-axis observations, the angular resolution is 50 arcseconds (50% 
encircled-photon radius) , and the maximum on-axis effective area of the system is about 260 cm^ 
(Pfeffermann et al. 1987). 

Spectra were constructed using the PSPC data for four spatially separate regions: the luminous 
northwestern and southwestern rims, the northeastern rim and the interior region. This final region 
was considered in order to analyze diffuse emission from the interior of the SNR. Emission from the 
two bright interior point sources - IWGA J1714.4-3945 and IWGA J1713.4-3949 (Pfeffermann 
&; Aschenbach 1996; Slane et al. 1999) - were subtracted from this region before a spectrum was 
extracted. The former source appears to be associated with a foreground star while the latter may 
be the X-ray counterpart to a recently-discovered radio pulsar. This radio pulsar may in turn be 
associated with G347.3— 0.5: we will discuss this source in more detail in Section 4.4. An image 
from the PSPC observation appears in Figure 1 with the positions of these regions indicated. Fits 
to the spectra extracted from these regions will be described in Section 3.2. 

2.2. ASCA GIS Observations and Data Reduction 

The instruments on ASCA (Tanaka, Inoue &: Holt 1994) observed locations on the southwestern 
rim, the northwestern rim and the northeastern rim of G347.3— 0.5. While both the GIS and the 
Solid-State Imaging Spectrometer (SIS) aboard ASCA collected data during these observations, 
the GIS has a larger field of view so only the data from the GIS were considered for the present 
study. The GIS consists of two units (denoted as GIS2 and GIS3) of imaging gas scintillation 
proportional counters with a sealed-off gas cell equipped with an imaging phototube: for our 
research, we only considered data collected by the GIS2 unit. The GIS is sensitive over the energy 
range of approximately 0.7 - 10 keV and has an energy resolution of AE/E = 0.078 at 6 keV 
(Ohashi et al. 1996; Makishima et al. 1996). The field of view of the GIS is 22' in radius and the 
maximum effective area for our off-axis observations is 94 cm^ at 2.07 keV. The data from these 
observations were reduced using standard ASCA GIS reduction techniques as described in "The 
ASCA Data Reduction Guide, Version 2.0." The reduction process filtered data for elevation angle, 
stable pointing directions, the South Atlantic Anomaly and cut-off rigidity. The data were also 
screened for characteristics of the GIS internal background and screening based on event locations 
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and rise time. The standard script "ASCASCREEN" was run in order to accomplish much of this 

reduction. We extracted spectra for the northwestern rim, the southwestern rim and the northeast 
rim of the SNR using the same regions as in the ROSAT PSPC spectral analysis. Fits to these 
spectra will be discussed in detail in Section 3.2. 

2.3. RXTE PCA Observations and Data Reduction 

G347.3— 0.5 was observed by the RXTE PCA with a pointing position that was offset by 
7.8' from the nominal center of the SNR to avoid contaminating flux from three nearby X-ray bina- 
ries (4U1708-408, IRXS J170849.0-400910 and GPS1713-388) that appear close to G347.3-0.5 
in the sky. Each of these X-ray binaries is located approximately 1° from the pointing position, 
so none of them appeared within the field of view (1° FWHM) of the RXTE PCA observation. 
Therefore, effects of contaminating flux from these X-ray binaries may be safely ignored. 

The PCA is a spectrophotometer comprised of an array of five coaligned proportional counter 
units that are mechanically collimated to have a field of view of 1° FWHM. The array is sensitive 
to photons that have energies between approximately 2 and 60 keV. The energy resolution AE/E 
of the array is 0.18 at 6 keV, and the maximum on-axis collecting area is about 6000 cm^ at 9.72 
keV. The PCA data were screened to remove the time intervals during which (1) one or more of 
the five proportional counter units is off, (2) G347.3— 0.5 is less than 10° above the limb of the 
Earth, (3) the background model is not well-defined and (4) the pointing direction of the detectors 
is more that 0°.02 from the nominal pointing direction in either right ascension or declination. 
After screening the data using these criteria, a final dataset was prepared. The PCA background 
spectrum for G347.3— 0.5 was estimated using the the FTOOL pcabackest^. 

G347.3— 0.5 lies toward the direction of the Galactic Ridge (GR), a region of the Galaxy which 
features extensive diffuse X-ray emission as described by Valinia &; Marshall (1998). G347.3— 0.5 
falls within the sub-region of the GR denoted as Rl by those authors: for this sub-region, Valinia 
&; Marshall (1998) obtained a best-fit for the observed X-ray emission using a two component 
model. The first component is a Raymond-Smith plasma component (Raymond & Smith 1977) - 
denoted as RAYMOND - which describes the emission spectrum from a hot diffuse gas, including 
line emissions from several elements. Using this model, with abundances frozen to solar, Valinia 
&; Marshall (1998) obtained a temperature of kT = 2.9 keV for the thermal emission from this 
sub-region of the GR. The second component was a simple power law (denoted as POWER LAW) 
which models the non-thermal emission from this sub-region: Valinia & Marshall (1998) obtained a 
photon index of 1.8 for this component. Lastly, the derived column density for the spectral fit was 
Nh = 1.8x10^^ cm^^. To sample this background X-ray emission, two additional observations were 
made with the PCA of fields adjacent to the field observed in the pointed observation of G347.3— 0.5: 



*See http://heasarc.gsfc.nasa.gov/docs/software/ftools/ftools_menu.html for more information about FTOOLS. 
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we denote these fields as Background Region 1 and Background Region 2, respectively (see Table 
1). We note that the approximate pointing errors for all of the PCA observations are ~8". Data 
from these two observations were reduced in the same way as the data from the pointed observation 
of G347.3— 0.5. In Figure 2, we present spectra extracted for the on-source observation and both 
of the background observations. We will return to a discussion of X-ray emission from the GR in 
Section 4.3. 



3. Spectral Fitting and Analysis 

3.1. The RXTE Spectrum of G347.3-0.5 

After reducing all of the data, we proceeded to perform spectral-fitting analysis using the X-ray 
spectral- fitting package XSPEC^, first considering the RXTE PCA data. Given that the X-ray 
emission from G347.3— 0.5 is known to be almost entirely non-thermal and probably produced by 
electrons emitting synchrotron radiation, we used three models for non-thermal X-ray emission to 
fit the data: a simple power law, the SRESC model and the SRCUT model. An emission feature 
is seen in the RXTE PCA spectrum near 6.4 keV: to model this feature, we included a Gaussian 
component with each of the models used in the spectral fitting (see Section 4.3 for a more detailed 
discussion about this apparent feature). In all cases, photoelectric absorption along the line of 
sight was modeled using the Wisconsin (Morrison &; McCammon 1983) cross-sections and relative 
abundances of elements as described by Anders & Ebihara (1982). 

POWER LAW: Earlier studies of the X-ray spectra associated with the bright rims of G347.3— 0.5 
by Koyama et al. (1997) and Slane et al. (1999) revealed that the X-ray emission from the north- 
western rim, the southwestern rim and the eastern region of G347.3— 0.5 could all be adequately 
fit using a simple power law (with a photon index of ~ 2). These fits were made based on ROSAT 
and ASCA observations (extending to approximately 10 keV), and we performed a fit to the RXTE 
data using a power law with approximately the same index. We included a second component, a 
Gaussian with a line width a frozen to keV, to fit the apparent emission feature seen near 6.4 
keV. We discuss this emission feature in more detail in Section 4.3. 

Figure 3 shows a POWER LAW-I-GAUSSIAN fit to the spectrum extracted from the RXTE 
PCA observations, using an photon index of 2.33 to correspond to the weighted mean value of pho- 
ton indices obtained in the fits made to the northwestern rim, southwestern rim and eastern region 
of G347.3— 0.5 by Slane et al. (1999). Clearly, this model did not fit the data in this energy range, 
yielding a reduced of 482. Even when the photon index was a free parameter, an acceptable fit 
was still not obtained: in Figure 4, we plot the results of a POWER LAW-hGAUSSIAN fit where 
we fit the photon index. In this Figure, the photon index is 2.61, yielding a of 1808 for 49 degrees 



®See http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/index.html for more information about XSPEC. 
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of freedom (a reduced of 36.9). While this is certainly an improvement, this model still does 
not fit the data. We conclude that the power law models employed by Slane et al. (1999) over the 
energy ranges sampled by ROSAT and ASCA cannot be extended into the RXTE energy range. 
The energy distribution of the electrons that produce the observed X-ray spectrum is therefore 
inconsistent with a simple power law. 

SRESC: The second model used in our analysis was SRESC, which describes a synchrotron 
spectrum from an electron distribution limited by particle escape above some energy. This model 
was used by Dyer et al. (2001) to sucessfuUy fit the X-ray spectrum of SN 1006, and thorough 
descriptions of this model are presented by Reynolds (1996) and Reynolds (1998). The SRESC 
model describes electrons which are shock-acclerated in a Sedov blast wave encountering a constant- 
density medium containing a uniform magnetic field. This model also includes variations in electron 
acceleration efficiency with shock obliquity and post-shock radiative and adiabatic losses. 

In Figure 5, we present a fit to the RXTE spectrum using the SRESC model combined with the 
GA USSIAN model. Similar to our attempts to fit the spectrum with a POWER LAW+ GA U SSI AN 
model, we could not obtain an acceptable fit to the RXTE spectrum of G347.3— 0.5 using this model. 
For the plotted model, the corresponding ratio of to degrees of freedom is 2451/48 = 51.07. We 
therfore conclude that the SRESC model, like the POWER LAW model, does not properly describe 
the non-thermal X-ray emission seen from this SNR. 

SRCUT: The third model we considered is denoted as SRCUT: it was described by Reynolds 
(1998), Reynolds & Keohane (1999) and Hendrick &; Reynolds (2001), with the latter two papers 
implementing the model in their study of the maximum energies of electrons accelerated by samples 
of SNRs in the Galaxy and the Large Magellanic Cloud, respectively. The SRCUT model describes 
a synchrotron spectrum from an exponentially cut off power law distribution of electrons in a 
uniform magnetic field. The photon spectrum is itself a cut-off power law, rolling off more slowly 
than an exponential in photon energies. Though this model is an oversimplification, it is more 
realistic than a power law and it does give the maximally curved physically plausible spectrum. 
This spectrum can in turn be used to set limits on maximum accelerated electron energies even 
in remnants whose X-rays are thermal. The SRCUT model assumes an electron energy spectrum 
A^e(^) of the form 

NeiE) = K E-^ e~^/^-^^°ff (1) 

where K'lsa, normalization constant derived from the observed fiux density of the region of the SNR 
at 1 GHz, r is defined as 2q;-|-1 (where a is the radio spectral index) and E cutoff is the maximum 
energy of the accelerated cosmic-ray electrons. A crucial advantage of this model (as well as the 
SRESC model) is that a resulting fit may be compared with two observable properties of an SNR, 
namely its flux density at I GHz and a. Finally, one of the fit parameters for the SRCUT and 
SRESC models is the cutoff frequency Vcutoff of the synchrotron spectrum of the electrons. This 
quantity is defined as the frequency at which the flux has dropped by a factor of 10 from a straight 
power law. We can express Vcutoff in a way which is quantitatively consistent with previous work 
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(Reynolds 1998; Reynolds & Keohane 1999; Hendrick & Reynolds 2001) as 

0.48 X 10- H. ,2) 

where B^g is the magnetic field strength of the SNR in /xG, assuming the electrons are moving 
perpendicular to the magnetic field. By using the value for I'cutoff returned by SRCUT as well 
as the normalization factor the synchrotron spectrum of the shock-accelerated electrons can 
be adequately described. Moreover, an estimate for the maximum energy E^utoff for the shock- 
accelerated electrons can also be derived. 

In Figure 6, we present a fit to the RXTE spectrum using the SRCUT model. Clearly, this 
model gives the best-fit to the spectrum of G347.3— 0.5 over this energy range: the corresponding 
value for x^/degrees of freedom is 207.98/48 = 4.33, which is a considerable improvement (though 
still not statistically acceptable) over the fits obtained by both the POWER LAW model and the 
SRESC model. The normalization K for the SRCUT component of this fit indicates that the 
flux density at 1 GHz for the entire SNR is 6.9 ± 0.2 Jy (90% confidence limits). This value is 
marginally consistent with the weak radio emission seen from this SNR, known to be dominated 
by a source of emission at the northwestern rim of G347.3— 0.5 with a corresponding flux density 
of 4 zb 1 Jy (Ellison et al. 2001). Moreover, the broad field of view of the PCA certainly intercepts 
diffuse X-ray flux from other sources besides G347.3— 0.5, which leads to a corresponding increase 
in the value for K returned by the flt. We also note that by using the best fit value for i'cutoff, 
Equation 2 yields E cutoff ~ 68 TeV (assuming a magnetic field strength of 10 ^G). This is in good 
agreement with the value of E cutoff obtained by Ellison ct al. (2001) using a similar value for B^^q 
with different models. Based on our success in fitting the RXTE spectrum of G347.3— 0.5 with the 
SRCUT model, we conclude that the energy distribution of the highest energy cosmic-ray electrons 
accelerated by this SNR is consistent with an exponentially cut-off power law. 

3.2. Joint Spectral Fitting of ROSAT PSPC, ASCA GIS and RXTE PCA Data 

We next applied the SRCUT model to the X-ray spectrum of G347.3— 0.5 as observed by the 
ROSAT PSPC and the ASCA GIS as weU as the RXTE PCA. We took advantage of the spatial 
resolving capabilities of the ROSAT PSPC and the ASCA GIS in order to simultaneously fit spectra 
extracted from particular regions of the SNR. We simultaneously fit eight different X-ray spectra: 
the spectra for the northwestern rim, the southwestern rim, the northeastern rim and the diffuse 
central emission (with the two interior point sources omitted) as observed by the ROSAT PSPC, 
the spectra for the northwestern rim, the southwestern rim and the northeastern rim as observed 
by the ASCA GIS and finally the spectra for the whole SNR as observed by the RXTE PCA. 

Two models (each a combination of a thermal component and a non-thermal component) 
were used for the joint spectral-fitting process - namely a RAYMOND+SRCUT model and an 
EQUIL+SRCUT model - in order to fit both the thermal and non-thermal components respec- 
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tively in the individual specta. In their study of the X-ray spectrum of G347.3— 0.5 based on 
ROSAT PSPC and ASCA GIS data, Slane et al. (1999) found no evidence for thermal emission 
from any portion of the SNR. We revisted this issue by testing whether the joint fit to the individual 
spectra would be improved by the inclusion of a thermal component. Both thermal models describe 
the emission spectrum from a hot diffuse gas, including line emission from several elements, that is 
in a state of collisional ionization equilibrium. Both of these models have a common set of param- 
eters: a plasma temperature kT, an abundance Z, a redshift z and an emission measure EM. For 
the purposes of the present work, we froze the values of Z and z to unity (i.e., solar abundance) 
and zero, respectively. The emission measure EM is defined as 



where d is the distance to the source in centimeters and Ug and nn are the electron and hydrogen 
densities, respectively, in units of cm~^. We can evaluate this expression as follows: we recall that 
we have assumed a distance of 6 kiloparsecs to G347.3— 0.5 and we note that the angular size of 
this SNR is 65' x 55' (Green 2001), therefore the corresponding linear radius is 1.65 x 10^° cm. 
Assuming a spherical volume for the SNR with a filling factor of one-quarter (that is, assumed 
only one-quarter of the total volume of the SNR is filled with thermal particles that produce the 
observed emission), as well as uniform values for ng and nn throughout the volume, then Equation 
3 yields 



We can estimate the ambient densities of hydrogen surrounding G347.3— 0.5 from the values ob- 
tained from the spectral fits for EM (assuming Ue = 1-2 uh)- We will continue this discussion 
when we describe our detection of thermal X-ray emission from G347.3— 0.5 in Section 4.2. 

We have modeled the background diffuse X-ray emission from the GR as observed by the 
RXTE PCA with the same two component model used by Valinia & Marshall (1998) mentioned 
earlier, namely a RAYMOND+POWER LAW model where the abundance Z is frozen to solar 
values and the redshift z is frozen to zero. We froze the values for kT, photon index and N h to 
be 2.9 keV, 1.8 and 1.8 x 10^^ cm~^, respectively: these are the same best-fit values measured by 
Valinia & Marshall (1998). The normalizations for the RAYMOND and POWER LAW components 
were left as free parameters. 

In Figures 7 through 9, we present data, fits and residuals for the spectra extracted from 
the ROSAT PSPC regions (Figure 7), the ASCA GIS regions (Figure 8) and the whole SNR as 
observed by the RXTE PCA (Figure 9) using the RAYMOND+SRCUT model. Similarly, in 
Figures 10 through 12, we present data, fits and residuals for the same set of spectra using the 
EQUIL+SRCUT model. In Tables 2 and 3, we present the fit parameters for the depicted fits: the 
quality of the two fits is comparable for the two models (a x^/degrees of freedom of 808.25/471 
= 1.72 for the RAYMOND+SRCUT model compared to a /degrees of freedom of 843.76/472 
= 1.79 for the EQUIL+SRCUT model for 501 pulse height amplitude bins). We emphasize that. 




(3) 
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for both two-component models, the fits successfully reproduce the observed radio properties of 
G347.3— 0.5, namely the two estimates for the flux densities at 1 GHz for the northwestern rim are 
both consistent with the observed value of 4 ± 1 Jy (Ellison et al. 2001). Also, the estimates for the 
flux densities at that frequency for the other portions of the SNR are very modest, also consistent 
with observations. In the next section we discuss these fit parameters in more detail. 



4.1. Cutoff Energies of Electrons Shock- Accelerated by G347.3— 0.5 

The electron energy Eg associated with synchrotron photons radiated at energy Ex is given 
by the relation 



We can use this equation and the values for f cutoff derived by our fits to estimate the maximum 
energy E^utoff of the electrons accelerated by G347.3— 0.5. We took four values for ^cutoff yielded 
by our fits to the emission from the northwestern rim and the diffuse central emission using the 
SRCUT+RAYMOND model and the SRCUT+EQUIL model. Using Equation 5, Ex = hvcutoff 
yields an estimate for E cutoff (assuming a magnetic field of 10 //G), and we list our calculated values 
for E cutoff in Table 4. These values range from approximately 19 TeV to 25 TeV and are consistent 
with the values obtained by Ellison et al. (2001). We agree with the results and arguments of Slane 
et al. (1999) and Ellison et al. (2001) that G347.3— 0.5 is not accelerating cosmic-ray electrons to 
the "knee" energy of 3000 TeV. In Section 5, we describe a more rigorous attempt to estimate the 
maximum energy of cosmic-ray electrons accelerated by this SNR, where we simultaneously consider 
radio. X-ray and TeV emission from the northwestern rim. We will employ a more sophisticated 
fitting method where the ambient magnetic field strength is allowed to vary. 



Previous X-ray observations of G347.3— 0.5 by Koyama et al. (1997) and Slane et al. (1999) 
found no evidence for thermal emission from any portion of this SNR. Slane et al. (1999) placed 
upper limits on the amount of thermal emission as a function of kT for both emission from the 
northwestern rim and the central diffuse emission of G347.3— 0.5. By including a thermal component 
in fits to the X-ray emission from G347.3— 0.5, we re-examine the upper limits on thermal emission 
derived by Slane et al. (1999). 

In Table 4, we present values for E cutoff-, T, uh and ng = 1.2 uh, as calculated from our model 

fit values for a, EM and Equation 4. We have considered only the fit vahics for the northwestern 
rim and the central diffuse emission based on the results of both the SRCUT+RAYMOND model 



4. 



Results 
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4.2. Thermal X-ray Emission from G347.3-0.5? 
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and the SRCUT+EQUIL model. For the northwestern rim, we can only provide upper limits 
for riH and n^- In Figures 13 and 14, wc plot our values for EM for the northwestern rim and 
the central diffuse emission, respectively, against the upper limits derived by Slanc ct al. (1999). 
We find that for the northwestern rim, our values for EM fall well below the upper limits of 
Slane et al. (1999), but for the diffuse emission form the central region, our values for EM more 
closely straddle the upper limits, and in the case of SRCUT+ RAYMOND our value exceeds the 
upper limit. This result indicates the first detection of thermal X-ray emission from G347.3— 0.5: 
we speculate that earlier work may have missed this thermal emission because earlier efforts had 
concentrated on the non-thermal properties of this SNR (particularly the bright rims of this SNR) 
and only considered a narrower energy range of the spectrum of G347.3— 0.5. In addition, the rather 
large angular size of this SNR makes it difficult to completely sample all of its X-ray properties. 
A more thorough understanding of the rather modest thermal emission from this SNR will require 
additional observations and modeling. 



4.3. The Detection of Iron Emission — Galactic Ridge Background? 

We now comment on a modest spectral feature seen near 6.4 kcV in the RXTE PCA spectrum 
of G347.3— 0.5 and at a similar energy in the spectra from the two background observations also 
made with the RXTE PCA. This emission feature may be associated with iron, which is known to 
produce a broad set of emission features between the energies of 6 keV and 7 keV. One possibility 
to consider is that the feature is produced by thermal emission from G347.3— 0.5: we note that 
Bykov (2002) has recently argued that fast moving isolated fragments of SN ejecta composed of 
heavy elements (such as iron) should produce Ka fiuoresence emission lines at X-ray energies, and 
that the iron Ka fiuoresence emission line energy is known to be approximately 6.4 keV. A second 
possibility is that the feature originates from diffuse emission associated with the GR, which is seen 
in projection beyond G347.3— 0.5. The large field of view of the RXTE PCA (approximately 1° 
FWHM) and the lack of imaging capabilities of this instrument make it difficult to clearly associate 
this feature with either G347.3— 0.5 or the GR. Recently, Tanaka (2002) presented an analysis 
of spectral features - including the iron Ka fiuoresence emission line - associated with different 
portions of the GR, as observed by the ASCA SIS. 

In Table 5, we present the line energies and normalizations (along with the 90% confidence 
limits for these values) for the emission feature as seen in the RXTE PCA spectra of G347.3— 0.5 
and the two background pointings. We also performed a simultaneous fit of the emission feature 
as seen in the three spectra, where the line energies and the normalizations for the feature as seen 
in each spectra were fixed to be the same. We also present the best-fit values for the line energy 
and the normalization for the joint fit as well. Notice that in each case, the line energies and the 
normalizations are all approximately the same, within the confidence limits. In order to check 
for any instrumental effects that may be corrupting the data, we obtained archived data from an 
RXTE PCA observation of Gas A made on 1999 August 5 (OBSID 40806-01-04-00): we selected 
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this particular observation of Cas A for analysis because it took place nearly contemporaneously to 
our G347.3— 0.5 observation (less than two months later). Using the same data reduction techniques 
described in Section 2.3, we generated a source spectrum of Cas A and measured both the line energy 
and the normalization of the iron Ka emission line seen in the spectrum of Cas A. We measured 
a line energy of 6.66^^:^^ keV, which is in good agreement with the energy of this emission line in 
the ASCA SIS spectra of Cas A as observed by Holt et al. (1994). Thus, we conclude that there 
are no instrumental effects that are corrupting our analysis of the emission feature near 6.4 keV. 

We therefore hypothesize that the spectral feature seen in the RXTE PCA is not associated 
with G347.3— 0.5 for two major reasons. First, this feature was not detected by either Slane et 
al. (1999) or by us during the analysis of ASCA CIS observations of this SNR which sampled the 
energy range which includes this spectral feature. Second, by inspection of Figure 2 and Table 5, 
it is clear that the emission feature is seen in both the on-source RXTE PCA pointing and the 
two background pointings, and the line energies and normalizations of the emission features for all 
three pointings are similar. If this feature was associated with G347.3— 0.5 rather than the GR, we 
would expect that the line energy and/or the normalization in the source region would significantly 
differ from the values seen for the feature in the two background pointings. In the study of spectral 
features observed from diffuse X-ray emission from the GR, Tanaka (2002) found that only diffuse 
emission from within 1° of the Galactic Center exhibited a feature near 6.4 keV, while the diffuse 
line spectrum from two other regions - centered at longitudes of 10° and 28.5°, respectively - 
were dominated instead by a feature near 6.7 keV. Because neither of these two regions include 
the positions of G347.3— 0.5 or the two background regions, however, we argue that the results 
from Tanaka (2002) are not applicable to our study presented here. We therefore conclude that 
the spectral feature is not associated with G347.3— 0.5 and instead originates in the diffuse X-ray 
emission from the GR seen in projection beyond G347.3— 0.5. 

4.4. X-ray Pulsations from the Radio Pulsar PSR J1713-3949 (Possibly Associated 

with IWGA J1713.4-3949 and G347.3-0.5)? 

Recent radio observations of IWGA J1713.4— 3949 (the X-ray source centrally located inside 
G347.3-0.5) by Crawford et al. (2002) discovered a new pulsar (PSR J1713-3949) with a period 

of 392 ms within 7 arcminutes (the discovery beam radius) of the X-ray source. Crawford et 
al. (2002) derived a distance of 5.0±0.2 kiloparsecs to this source, which is consistent with the 
distance estimate to G347.3-0.5 of Slane et al. (1999). The position of IWGA J1713.4-3949 near 
the center of G347.3— 0.5 suggests that the X-ray source (as well as the radio pulsar) and the SNR 
are associated. In addition, the position of G347.3— 0.5 itself near a molecular cloud also suggests 
that the SNR was produced by a massive star in a core collapse SNe: these types of SNe are 
expected to produce neutron stars. To further investigate the nature of IWGA J1713.4— 3949, we 
performed a timing analysis of our RXTE PCA data to search for pulsations from this source. 

In Figure 15, we present a power spectrum of our RXTE PCA data over the frequency range 
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of 2.5 through 2.6 Hz with a dashed vertical hne indicating a period of 392 ms (2.55 Hz). We find 

no evidence for pulsations in our spectrum at any frequency, including 2.55 Hz. Our failure to 
detect pulsations from this central source may be explained by several factors, such as significant 
contaminating diffuse background emission from the GR, cosmic ray background and finally emis- 
sion from the SNR dominating over emission from a pulsar over the energy range sampled by our 
RXTE PCA observation. Pointed high angular resolution X-ray observations would be much more 
sensitive to pulsed X-ray emission from this source. 



5. The Broadband Energy Spectrum of the Northwestern Rim of G347.3— 0.5 

In Figure 16, we present a broadband energy spectrum of the northwestern rim of G347.3— 0.5 
ranging from radio through 7-ray energies. We have constructed this diagram using published values 
for the 1 GHz radio observation of G347.3— 0.5 by the Australian Telescope Compact Array (ATCA) 
as presented by Ellison et al. (2001), the ROSAT PSPC and ASCA GIS data described previously 
in this paper and data from the CANGAROO observation at TeV energies of the spectrum of 
this rim as published by Enomoto et al. (2002). Using the ISIS software package^ (Houck & 
Denicola 2000), we obtained a fit to the energy spectrum using a model that includes multiple high 
energy processes associated with SNRs, as indicated by separate curves. These curves indicate 
the emission profile expected by a shock-accclcratcd population of electrons from two different 
processes: synchrotron emission (as indicated by the curve labeled "S"), and inverse-Compton 
emission from cosmic microwave background photons that have been upscattered to TeV energies 
by energetic cosmic-ray electrons (the curve labeled "IC"). These emission models are consistent 
with the models employed by Sturner et al. (1997), and a detailed discussion of the high energy 
emission from non-linear shock acceleration processes associated with SNRs is given by Baring et 
al. (1999). From our fit to the ASCA GIS data and the TeV data, we have derived estimates for 
E cutoff for the accelerated cosmic-ray electrons and magnetic field strength B. For the depicted fit, 
the ratio of to the degrees of freedom is 191.78/113 = 1.70. Based on this fit, we estimate a value 
of only 8.8I34 TeV for E cutoff of the accelerated cosmic-ray electrons and a corresponding magnetic 
field strength oi B = 1501|q^ /iG. This value for E cutoff is sharply lower (but still consistent) with 
the upper limit derived by Ellison et al. (2001), while the corresponding value for B is considerably 
larger than the value of 10//G assumed by Slane et al. (1999). There is also a modest amount of 
overlap between our range of values for B and the the range suggested by Enomoto et al. (2002) (10 
/xG - 100 /xG). However, the physical plausibility of this fit may be questioned based on the value 
derived for the ratio of the volumes for TeV emission and X-ray emission from the northwestern rim 
— that is, VtcV /V X-ray In contrast to Enomoto et al. (2002), who fixed this ratio to equal unity, 
we allowed this ratio to vary and obtained a value of VtcV x-ray ~ 1000, with a lower limit of 
360, assuming that the power law index is 2.1 (derived from our from our model fits where a = 



See http://space.mit.edu/ASC/ISlC/ for more information about ISIS. 
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0.55) and that there is no curvature to the electron spectrum. This ratio is too high to be physically 

acceptable, and suggests that the TeV emission from the northwestern rim of G347.3— 0.5 is not 
produced by inverse Compton scattering of cosmic microwave background photons off high energy 
electrons accelerated along this rim of the SNR. 

Nonetheless, while this TeV emission may not be produced by inverse Compton scattering, we 
argue that neutral pion decay also cannot produce this emission. We point out that in arguing that 
the TeV emission is produced by neutral pion decay, Enomoto et al. (2002) assumed an extremely 
large ambient density (n pa 100 cm~^ if a distance of 6 kpc to G347.3— 0.5 is assumed) for the 
northwestern rim of the SNR. Our fits to the thermal component of the X-ray emission from the 
northwestern rim and the central diffuse emission yield a dramatically lower ambient density (n 
~ 0.05-0.07 cm~^) and make the conclusions of Enomoto et al. (2002) seem untenable. Slane 
et al. (1999) also concluded that G347.3— 0.5 is expanding into an ambient interstellar medium 
with an extremely low density, and suggested that the SNR is expanding into the evacuated wind- 
blown cavity created by the stellar progenitor: we agree with this finding, and suggest that while 
G347.3— 0.5 may have formed in the vicinity of dense molecular clouds, it has not yet begun to 
interact with these clouds. We also note that the observed TeV emission is observed to be centered 
on the northwestern rim of the SNR and not on the adjacent cloud, as assumed by Enomoto et al. 
(2002). Finally, Reimer &: Pohl (2002) noted that if the flux of the GeV emission from the 7-ray 
source 3EG J1714— 3857 can be used to constrain the neutral pion decay spectrum fit presented 
by Enomoto et al. (2002): if this source is associated with the northwestern rim of G347.3— 0.5, 
its flux mTist be taken into account in all spectral fits, while if the 7-ray source is not associated 
with this rim, then the GeV emission may be treated instead as upper limits on flux along this 
energy range. Reimer &: Pohl (2002) showed that the flux from the neutral pion decay spectrum 
fit exceeded the observed GeV emission by a factor of three, therefore strongly arguing against 
neutral pion decay as the process responsible for the TeV emission. For these reasons, we argue 
that the TeV emission observed from the northwestern rim of G347.3— 0.5 is not produced by 
neutral pion decay. We also comment that if the ambient density is as low as inferred by our 
estimates, the process of non-thermal bremsstrahlung emission - as well as neutral pion decay - 
will be significantly suppressed. A non-thermal bremsstrahlung origin for the TeV emission may 
also be ruled out because a spectrum generated by this process would - like the spectrum produced 
by neutral pion decay - exceed the observed GeV emission. In conclusion, we state that none 
of the known high energy processes associated with SNRs - inverse Compton scattering, neutral 
pion decay and non-thermal bremsstrahlung - can adequately fit the observed TeV emission from 
the northwestern rim. From the work of Ellison et al. (2000), we suspect that all models of high 
energy emission from SNRs over broad energy ranges must take into account some curvature in 
the energy spectrum of accelerated electrons rather than assuming simple power laws. In such a 
case, one of these emission processes may be able to adequately fit the observed TeV emission from 
the northwestern rim of G347.3— 0.5 - and perhaps the rims of other SNRs detected at such high 
energies - after all. We will discuss such curvature in the energy spectra of SNRs in more detail in 
a future work. 
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Lastly, we consider the putative association between G347.3— 0.5 and the X-ray source AX 
J1714. 1—3912, a hard spectrum source coincident with a molecular cloud located along the north- 
eastern rim of G347.3— 0.5. Uchiyama et al. (2003) presented an analysis of the spectral properties 
of this source using data from ASCA GIS observations. The X-ray spectrum of this source was best 
fit by a power law with a rather fiat photon index of r=1.0lQ3, prompting Uchiyama et al. (2003) 
to argue that such a spectrum was best interpreted as non-thermal bremsstrahlung emission from 
particles accelerated at the shock interface between the northeastern rim of G347.3— 0.5 and the 
molecular cloud. Despite the spatial proximity between G347.3— 0.5 and the molecular cloud (and, 
by extension, AX J1714. 1—3912), the evidence for the association between these sources is not 
compelling because, based on simulations of the interactions between SNRs and interstellar clouds 
(Jun & Jones 1999), we expect the amount of radio emission from both the SNR and the cloud to 
dramatically increase as a result of the interaction. Jun k. Jones (1999) also predict that if the SNR 
has overtaken the cloud, the radio emission from the SNR would "wrap around" the cloud, and such 
a morphology should be readily apparent. Radio images of G347.3— 0.5 and its surrounding envi- 
ronment have been made with the Molonglo Observatory Synthesis Telescope (MOST) (Roberts 
1994) and presented by Slane et al. (1999) (at the frequency of 843 MHz) and Ellison et al. (2001) 
(at the frequency of 1.4 GHz). Inspection of these images reveals little (if any) radio emission at 
either the position of AX J1714. 1—3912 or the northeastern rim of G347.3— 0.5, and no evidence 
for a "wrap around" morphology in the radio emission. The failure to detect radio emission from 
the hard X-ray source or the northeastern rim of the SNR and a "wrap around" morphology to the 
radio emission argues against an interaction taking place at that site. Another reason to doubt that 
an interaction is taking place between G347.3— 0.5 and the molecular cloud is related to the ratio 
of emission from rotational transitions of CO molecules, C0(J=2^1)/C0(J=1— >^0), as measured 
at the site of the cloud. An enhanced value for this ratio (such as a value of 1.2 or greater) is 
thought to indicate that an interaction is taking place between the cloud and a shock, such as the 
one associated with an SNR (Seta et al. 1998). Both Slane et al. (1999) and Butt et al. (2001) 
found an enhanced value for this ratio for the cloud along the northeastern rim of the SNR, but 
we suggest that in this particular case the enhanced CO(J=2^1)/CO(J=1^0) ratio and the low 
ambient density of the SNR cannot be reconciled, and that (as mentioned before) G347.3— 0.5 is 
not currently interacting with this molecular cloud. Butt et al. (2001) measured a value for the 
flux ratio CO(J=2^1)/CO(J=1^0) of 2.4±0.9: according to Draine & Roberge (1984), such a 
flux ratio corresponds to a shock velocity of less than 5 km/sec in the cloud. The shock speed in a 
cloud is reduced by the square root of the ratio of the ambient densities. Based on CO observations 
(Bronfman et al. 1989; Butt et al. 2001), the mean density of the cloud has been estimated to 
be approximately 500 cm~^. Prom our spectral fits, we have inferred an average ambient density 
G347.3— 0.5 to be ~ 0.06 cm~^: from these estimates for the mean density of the cloud and ambient 
density around the SNR, we conclude that the shock velocity outside the cloud is less than 500 
km/sec. For comparison, if we assume that the electrons and protons in G347.3— 0.5 are in thermal 
equilibrium at kT = 1.5 keV (from our fits to the diffuse thermal emission from the SNR), then 
the corresponding shock velocity is approximately 1100 km/sec. Since the proton temperature is 
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most likely much larger than the electron temperature in the case of young SNRs (Hwang et al. 
2002), like G347.3— 0.5, this value for the shock velocity is actually a lower limit, and therefore the 
shock velocity exceeds 1100 km/sec, which is in sharp disagreement with the shock velocity of < 
500 km/sec derived earlier. Uchiyama et al. (2003) noted that the large X-ray luminosity of AX 
J1714.1-3912 (1.7 X 10^^ ergs/sec at the assumed distance of 6 kpc) made G347.3-0.5 barely able 
to provide enough kinetic energy (in the form of nonthermal particles) to power this source, unless 
G347.3— 0.5 is much closer than 6 kpc or the SN explosion which created this SNR had more than 
Iq51 gj,gg q£ energy. We suggest that G347.3— 0.5 is located in the vicinity of the clouds, but has 
not yet started to interact with any of the clouds. More observations and analysis are required to 
determine the true nature of AX J1714. 1—3912. 

6. Comparison with Other Shell- Type SNRs with Non-Thermal X-ray Emission: 

SN 1006 and G266.2-1.2 

In Table 6, we compare the gross properties of G347.3— 0.5 with two other dynamically-young 
shell-type SNRs which possess significant non-thermal components to their X-ray emission, SN 1006 
(Dyer et al. 2001; Allen, Pctrc & Gotthelf 2001) and G266.2-1.2 (Slane et al. 2001). Inspection 
of this Table 6 reveals that all three SNRs are X-ray luminous, radio faint and are expanding into 
regions of the interstellar medium with particularly low ambient density. Both core-collapse and 
Type la SNe appear to produce SNRs with non-thermal X-ray emission: SN 1006 is the archetype 
of a SNR produced by a Type la SN, while G266.2-1.2 and G347.3-0.5 are thought to be produced 
by Type II SNe. The X-ray morphologies of these three SNRs also contrast as well: for SN 1006, 
the morphology is bilateral with two X-ray bright rims, visible in ASCA images of this SNR (Dyer 
et al. 2001), while the morphologies of both G266.2— 1.2 and G347.3— 0.5 feature multiple bright 
rims: ROSAT and ASCA images of these two SNRs reveal luminous northwestern, southwestern 
and northeastern rims for both SNRs, while G266.2— 1.2 also has a luminous southern rim as well 
(Slane et al. 1999, 2001). The X-ray synchrotron flux depends on the exponential cut-off energy, 
the strength of the magnetic field, and the density of nonthermal electrons. In the case of SN 1006, 
the morphology may be due to an enhancement (via compression) in the magnetic field along the 
two bright rims relative to the magnetic field strength elsewhere along the rim. The morphology of 
SN 1006 - with two symmetric luminous X-ray rims - was reproduced remarkably well by the work 
of Reynolds (1998), who described models of synchrotron X-rays from shell supernova remnants. In 
the cases of G347.3— 0.5 and G266.2— 1.2, the cause of the morphology is not clear. The bright rims 
may be associated with regions where the density of nonthermal electrons is relatively high. The 
common property of low ambient density for these SNRs does not imply an exclusive requirement for 
non-thermal X-ray emission: for example, another SNR, Gas A, possesses a strong non-thermal X- 
ray component to its emission (Vink &; Laming 2003) - though not a component which dominates 
the X-ray spectrum of the SNR - and is expanding into a particularly dense ambient medium. 
Nonetheless, it appears that a low ambient density - independent of the type of SN progenitor - is 
conducive to the X-ray spectrum of the SNR being dominated by non-thermal X-ray emission. Deep 
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X-ray observations and detailed analysis of other SNRs located in regions of low ambient density 

are required to determine the degree to which low ambient density is necessary for producing this 
type of X-ray emission: we will explore this subject in more detail using a sample of such SNRs in 
a future work. 

7. Conclusions 

The results and conclusions of this work may be summarized as follows: 

1) We present a spatially-resolved X-ray spectral analysis of both the three bright rims and 
the diffuse central emission of the Galactic SNR G347.3— 0.5. This analysis involves data from 
observations made of this source using instruments aboard three different X-ray satellites (namely 
the ROSAT PSPC, the ASCA GIS and the RXTE PCA) spanning the approximate energy range 
of 0.5 - 30 keV. 

2) We have successfully fit the spatially-resolved X-ray spectra of G347.3— 0.5 using the SRCUT 
model, describing an exponential cut-off in the population of the highest energy electrons accelerated 
by this SNR. In contrast, neither the SRESC model nor a simple power law model yielded acceptable 
fits to the data. Wc find that the power law fits obtained by Slanc ct al. (1999) to fit the spectrum 
of G347.3— 0.5 over the energy range of 0.5-8 keV do not adequately fit the spectrum of this SNR at 
higher energies. We have also detected for the first time a thermal component to the X-ray emission 
from this SNR: this component appears to be more closely associated with central diffuse emission 
from the SNR than from the X-ray luminous rims. We estimate the ambient density surrounding 
the SNR to be f«0.05-0.07 cm""^: this value is consistent with the range of values obtained by Slane 
et al. (1999). 

3) A weak emission feature is seen near 6.4 keV in our RXTE PCA spectra of G347.3— 0.5 
and our two background pointings which sample the GR. Because the strength and the line energy 
of this feature are approximately the same in all three pointings, we argue that the feature is 
associated with diffuse background emission from the GR rather than from G347.3— 0.5 itself. 

4) We have analyzed our RXTE PCA data to search for X-ray pulsations from the radio 
pulsar PSR J1713— 3949, which lies near the center of G347.3— 0.5 and may be associated with the 
SNR. We cannot confirm the presence of X-ray pulsations from this source, though our data is 
dominated by both cosmic ray background and X-ray emission from the SNR and the GR, making 
the detection of pulsed signal from the pulsar difficult. 

5) Based on our two component (thermal and non-thermal) fits to the X-ray emission from 
the northwestern rim and the central diffuse emission from G347.3— 0.5 and assuming a magnetic 
field strength of i? = 10 fiG, we estimate the maximum energy E cutoff of accelerated cosmic ray 
electrons to be 19-25 TeV, consistent with previous analyses. Fitting the broadband (radio to 7- 
ray) energy spectrum of G347.3— 0.5 with a synchrotron-inverse Compton scattering model yields 
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values of 8.81^;^ TeV for the maximum energy E cutoff of the accelerated cosmic-ray electrons and 
a magnetic field of B = fiG. Our value for Ecutoff is sharply lower (but still consistent) 

with upper limits derived by previous studies, while our value for B is significantly higher than 
previously assumed. However, the ratio of volumes of TeV emission to X-ray emission derived by 
this fit {VTev/V X-ray ~ 1000, with a lower limit of 360) is too large to be physically reasonable. 
While it appears that inverse Compton scattering cannot adequately fit the TeV emission from 
the northwestern rim of G347.3— 0.5, the other two high energy processes associated with SNRs - 
neutral pion decay, as argued by Enomoto et al. (2002), and non-thermal bremsstrahlung - can 
also be ruled out as well. This situation may change with the application of more sophisticated 
models which take into account curvature of the energy spectrum of electrons. 

6) We have considered the gross properties of G347.3— 0.5 and two other SNRs known to 
feature X-ray spectra dominated by non-thermal emission, SN 1006 and G266.2— 1.2. We find that 
all three of these SNRs are dynamically young. X-ray luminous but radio weak SNRs expanding 
into regions of low ambient density. This suggests that low ambient density may play a key role 
in dictating that the X-ray spectrum from a SNR is dominated by non-thermal emission, though a 
more detailed study of the X-ray spectrum of other SNRs is required in order to test this hypothesis. 
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Fig. 1.— ROSAT PSPC image of G347.3-0.5 from Slane et al. (1999). Note the X-ray luminous 
northwestern rim (a TeV source) as well as the bright southwestern and northeastern rims. Contour 
levels start at 1.44 counts per arcmin"^ sec~^ and increase upward in steps of 1.15 counts per 
arcmin"^ sec"-*^. The eastern bright source (labeled as "#2") is a foreground star, while the central 
source (labeled as "#1" and also known as IWGA J1713.4— 3949) may be associated with the 
radio pulsar PSR J1713— 3949. Discrete regions of extracted spectra - namely the northwestern 
rim ("NW"), the southwestern rim ("SW") and the northeastern rim ("NE") - are indicated by 
solid boxes. The central diffuse emission ("C") is indicated by a dashed box: flux from the three 
rims and the two interior point sources were excised from the flux from this region. 
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Fig. 2. — Spectra from RXTE PCA observations of G347.3— 0.5 and two background regions (la- 
beled as "BGD#1" and "BGD#2") using an SRCUT model combined with a Power Law+Raymond 
component to fit the Galactic Ridge emission using model parameters derived by Valinia & Marshall 
(1998). 
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Fig. 3. — Power Law (Frozen Index) + Gaussian fit to the RXTE PCA observations of G347.3— 0.5. 
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Fig. 4 — Power Law (Free Index) + Gaussian fit to the RXTE PCA observations of G347.3-0.5. 



-27- 



100.00 



RXTE PCA -- SRESC+GAUSSIAN 



> 

CD 
C/) 

o 
O 



10.00 



1.00 



0.10 



^0.01 

^ 2.0 
^ 1.5 
B 1-0 

Q 




5 10 
Energy (keV) 



20 30 



Fig. 5.— SRESC + Gaussian fit to the RXTE PCA observations of G347.3-0.5. 
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Fig. 6.~ SRCUT + Gaussian fit to the RXTE PGA observations of G347.3-0.5. 
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Fig. 7. — Spectra of the center (C), the northwestern rim (NW), the northeastern rim (NE) and 
the southwestern rim (SW) of G347.3— 0.5 as observed by the ROSAT PSPC and fit using the 
RA YMOND+ SR CUT model. 
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Fig. 8. — Spectra of the northwestern rim (NW), the southwestern rim (SW) and the northeastern 
rim (NE) of G347.3-0.5 as observed by the ASCA GIS and fit using the RAYMOND+SRCUT 
model. 
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Fig. 9.— Spectrum of G347.3-0.5 as observed by the RXTE PGA and fit using the RAY- 
MOND+SRCUT model. 
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Fig. 10. — Spectra of the center (C), the northwestern rim (NW), the northeastern rim (NE) and 
the southwestern rim (SW) of G347.3-0.5 as observed by the ROSAT PSPC and fit using the 
EQUIL+SRCUT model. 
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Fig. 11. — Spectra of the northwestern rim (NW), the southwestern rim (SW) and the northeastern 
(NE) rim of G347.3-0.5 as observed by the ASCA GIS and fit using the EQUIL+SRCUT model. 
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Fig. 12. — Spectrum of G347.3— 0.5 as observed by the RXTE PCA and fit using the 
EQUIL+SRCUT model. 
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Fig. 13. — Upper limits on Emission Measure for the northwestern rim of G347.3— 0.5, as calculated 
by Slane et al. (1999). The values (corresponding to 90% confidence limits) for the emission measure 
from our fits to the northwestern rim using the RAYMOND+ SRCUT and EQUIL+SRCUT are 
plotted using a diamond and a square, respectively. 
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Fig. 14. — Estimates of the Emission Measure for the central diffuse emission from G347.3— 0.5, as 
calculated by Slane et al. (1999). The values (with 90% confidence intervals) for the emission mea- 
sure from our fits to the northwestern rim using the RAYMOND+ SRCUT and EQUIL+SRCUT 
are plotted using a diamond and a square, respectively. 
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Fig. 15. — Power spectrum of RXTEVCA. observation of G347.3— 0.5. The frequency corresponding 
to the putative pulsar (392 ms or 2.55 Hz) is indicated by the vertical dashed line. No evidence for 
pulsations is detected at this frequency. 
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Fig. 16. — Radio to gamma-ray photon energy-flux spectrum of the northwestern rim of 
G347.3-0.5. See Section 5. 
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Table 1. ROSAT PSPC, ASCA GIS and RXTE PGA Observations of G347.3-0.5 









Observed 






Exposure 


Satellite and 




Observation 


Portion of 


R.A. 


Decl. 


Time 


Instrument 


Obs ID 


Date 


G347.3-0.5 


{J2000.0) 


(J2000.0) 


(Seconds) 


ROSAT PSPC 


RP500198N00 


1992 September 23 


All 


17 13 33.6 


-39 48 36 


2758 


ASCA GIS 


55057000 


1997 March 23 


NW Rim 


17 12 17.8 


-39 35 31 


5428 


ASCA GIS 


55058000 


1997 March 24 


SW Rim 


17 12 53.8 


-39 54 23 


5776 


ASCA GIS 


55059000 


1997 March 25 


NE Rim 


17 14 28.3 


-39 35 26 


3924 


iJXT£; PGA 


40145-01 


1999 June 15 


All 


17 14 13 


-39 50 30 


38371 


RXTE PGA 


40145-02 


1999 June 14 


Background #1 


17 22 26 


-37 18 10 


8064 


iJXTE PGA 


40145-03 


1999 June 14 


Background #2 


17 05 24 


-42 20 32 


3229 
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Table 2. Best-Fit Parameters - SRCUT+RAYMOND ModeP 



Region 


Model 


Parameter (Units) 


Value 


90% Confidence 
Limits 


IN W 


WAnb 


IS H (.J-U cm j 


U. /4 


U. < i, 0. 1 1 




DJXLiU 1 


a. 


U.oo 


U.oi, u.oy 




DJxLiU 1 


^cutoff (10 tiz) 


i.4o 


i.oo, i.oy 




dKCU 1 


Normalization K (Jy) 


/I no 
4.Uz 


o.o9, 4.U9 




riA 1 M UiyJJ 


Ki {Ke\ ) 


i.io 






7? A VM^l AID 
TLA 1 ivl UI\JJ 


Emission Measure EM^ 


^ A KA \y 1 n— 3 
4.04 X iU 








(cm ) 






IN ill 


WAnD 


AT ('in22 „™— 2\ 

l\ H [w cm j 


U.Di 


U.4o, U. 70 




dKCU I 


a 


O.o7 


O.o4, O.OZ 




oKLj U 1 


'^cutoff (10 Hzj 


z.b ^ 


Z.4U, z.yo 




bKLU 1 


Normalization (Jy) 


4.07 X 10 


1 oo in — 2 

1.60 X iU , 










nil 




JxA 1 Ivl UI\JJ 


i,m ^^ro\T\ 


1 1 F^ 

i.io 


n on 1 vi 




JxA ilViUlyJJ 


Emission Measure EM^ 


D.ii X iU 


Q ncc V/ 1 n— 3 
o.Uo X iU , 






icm ) 




Q fiO V 1 n~3 
y.oy X lu 


bW 


WABb 


AT l-{Ci22 „™— 2\ 

IS H (10 cm j 


O.DO 


O.dO, O.do 




bKlyU 1 


a 


U.OU 


n /I c n c^Q 
U.4o, U.OO 




bKLU 1 


^cutoff (10 Hzj 


O /I Q 


oT f^n 
z.2r, z.oO 




bKLU 1 


Normalization (Jy) 


0. /Z 


(\ AO n no 
0.4o, U.yz 




itA 1 lVIUi\JJ 


KI \K<d\ ) 


o.oi 


CO /I nn 
z.oz, 4.UU 




KAi jyiUlyJJ 


Emission Measure EM^ 


Q V/ 1 n— 3 
O.DD X iU 


^ 1 n— 3 
Z.OO X iU , 






(cm j 




4.78 X 10"^ 


C 


WABS 


Nh (10^2 cm-2) 


0.49 


0.47, 0.51 




SRCUT 


a 


0.48 


0.47, 0.51 




SRCUT 


'^cutoff (10^^ Hz) 


2.41 


2.29, 2.52 




SRCUT 


Normalization (Jy) 


1.28 


0.98, 1.30 




RAYMOND^ 


kT (keV) 


1.56 


1.50, 1.61 




RAYMOND 


Emission Measure EM^ 


3.85 X 10-2 


3.65 X 10-2, 






(cm~^) 




4.15 X 10-2 


Background'^ 


WABS 


Nh (10^2 cm-2) 


1.8^^ 






POWER LAW 


Photon Index 


1.8^1 






POWER LAW 


Normalization 


2.72 X 10-3 


2.57 X 10-3, 
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Table 2 — Continued 



Region Model 


Parameter (Units) 


Value 


90% Confidence 
Limits 




(Photons keV~^ cm~^ 




2.84 X 10-3 




sec~^ at 1 keV) 






RAYMOND"" 


kT (keV) 






RAYMOND 


Emission Measure EM^ 


3.46 X 10-2 


3.38 X 10-2, 




(cm""^) 




3.55 X 10-2 



^X^/Degrees of Freedom = 808.25/471 = 1.72 using 501 PHA bins. 

Abundance frozen at 1.0 and redshift frozen at 0.0. 
'^See Section 3. 

'^Galactic diffuse background as observed by the RXTE PCA over the energy 
range of 2 through 30 keV. Nh, photon index and kT have all been frozen to the 
values measured by Valinia &; Marshall (1998). 



Table 3. Best-Fit Parameters - SRCUT+EQUIL Model^ 



90% Confidence 







JT dJ. dililc Lcl U 111 to J 


Value 


T 1 YYt tJ"C 


NW 


WABS 


Nh (10^2 cm-2) 


0.74 


0.70, 0.79 




SRCUT 


a 


0.55 


0.49, 0.57 




SRCUT 


i^cutoff (10^^ Hz) 


1.47 


1.31, 1.48 




SRCUT 


Normalization (Jy) 


4.48 


2.70, 4.55 




EQUIL^ 


kT (keV) 


0.49 






EQUIL 


Emission Measure EM'^ 
(cm~^) 


<4.64 X 10-3 




NE 


WABS 


Nh (10^2 cm-2) 


0.66 


0.56, 0.75 




SRCUT 


a 


0.44 


0.37, 0.50 




SRCUT 


cutoff (10^^ Hz) 


2.93 


2.59, 3.48 




SRCUT 


Normalization (Jy) 


0.14 


0.03, 0.51 




EQUIL^ 


kT (keV) 


1.34 


1.08, 1.71 




EQUIL 


Emission Measure EM'^ 
(cm~^) 


9.31 X 10-3 


6.56 X 10-3, 
1.31 X 10-2 


SW 


WABS 


Nh (10^2 cm-2) 


0.61 


0.58, 0.69 




SRCUT 


a 


0.52 


0.49, 0.55 




SRCUT 


'^cutoff (10^^ Hz) 


2.47 


2.31, 2.59 




SRCUT 


Normalization (Jy) 


1.10 


0.58, 1.74 




EQUIL^ 


kT (keV) 


1.31 






EQUIL 


Emission Measure EM'^ 
(cm~^) 


<2.52 X 10-3 




c 


WABS 


Nh (1022 cm-2) 


0.49 


0.48, 0.51 




SRCUT 


a 


0.50 


0.49, 0.52 




SRCUT 


cutoff (10^^ Hz) 


2.44 


2.34, 2.53 




SRCUT 


Normalization (Jy) 


1.94 


1.42, 2.61 




EQUIL^ 


kT (keV) 


1.69 


1.62, 1.74 




EQUIL 


Emission Measure EM'^ 

(cm"^) 


3.61 X 10-2 


3.45 X 10-2, 
3.82 X 10-2 


Background*^ 


WABS 


Nh (1022 cm-2) 


1.8<^ 






POWER LAW 


Photon Index 


1.8^ 






POWER LAW 


Normalization 
(Photons keV"^ cm-2 


3.07 X 10-3 


2.91 X 10-3, 
3.18 X 10-3 
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Table 3 — Continued 



Region Model 


Parameter (Units) 


Value 


90% Confidence 
Limits 




sec"-*^ at 1 keV) 






RAYMOND^ 


kT (keV) 


2.9^ 




RAYMOND 


Emission Measure EM'^ 


3.11 X 10-2 


3.04 X 10-2, 




(cm~^) 




3.21 X 10-2 



^X^/Degrees of Freedom = 843.76/472 = 1.79 using 501 PHA bins. 

^Abundance frozen at 1.0 and redshift frozen at 0.0. 
^See Section 3. 

'^Galactic diffuse background as observed by the RXTE PCA over the energy 
range of 2 through 30 keV. Nn, photon index and kT have all been frozen to the 
values measured by Valinia & Marshall (1998). 
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Table 4: Derived Values for E cutoffs T, nn and Uf, Based on Model Fit Parameters 

SRCUT + RAYMOND SRCUT + EQUIL 



Pariuuctor 




C 






E cutoff (TeVr 


19.6 


25.3 


19.7 


25.4 


pb 


2.10 


1.96 


2.10 


2.00 


uh (cm~^) 


<1.86 X 10-2 


0.054 


<1.89 X 10-2 


0.052 


Ue (cm~^) 


<2.23 X 10-2 


0.065 


<2.26 X 10-2 


0.062 



"Assuming a magnetic field of lO/uG. 
'T = 2a + 1. 
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Table 5. Line Energies and Normalizations for Emission Feature Detected Near 6.4 keV in 

RXTE PCA Spectra 





Line 


90% Confidence 




90% Confidence 


Pointing 


Energy (keV) 


Limits 


Normalization^ 


Limits^ 


G347.3-0.5 


6.41 


6.38, 6.42 


5.05 


4.70, 5.40 


Background ^1 


6.51 


6.47, 6.56 


4.80 


4.11, 5.59 


Background #2 


6.26 


6.13, 6.37 


7.24 


5.48, 9.91 


Joint Fit^ 


6.43 


6.40, 6.45 


5.06 


4.76, 5.38 


Gas 


6.66 


6.64, 6.67 


58.4 


57.8, 59.2 



^In units of 10 ^ photons cm ^ sec ^ . 

''Simultaneous fitting to the RXTE PCA spectra of G347.3-0.5 and both of the back- 
ground pointings. 

'^See Section 4.3. 
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Table 6. Properties of Shell- Type SNRs with Non-Thermal X-ray Emission 



I'ropcriy 


SN 1(K)() 


l^(lci(UiC(' 


G2()().2-1.2 


Hclcrcncc 


G:-!i7.:-!-0,5 




Distance (kpc) 


2.2 


(1) 


1-2 


(2) 


6 


(3), (4) 


Ambient Density N (cm~^) 


0.1 


(5) 


0.05a 


(2) 


?b0.05-0.07 


(6) 


Angular Size (arcmin) 


30 


(7) 


120 


(7) 


65 X 55 


(7) 


Dynamical Age (yr) 


996 


(1) 


ail4000=' 


(2) 


RiSOOO 


(4) 


Flux Density ^20 (Jy) 


19 


(7) 


50(?) 


(7) 


4 


(4) 


a 


0.6 


(7) 


0.3(?) 


(7) 


RiO.55 


(6) 


B (mG) 


40 


(5) 






150tl«° 


(6) 


B-field Direction 


Radial 


(8) 


Tangential 


(9) 






Progenitor Type 


la 


(1) 


II 


(2) 


II 


(3), (4) 



Note. — References: (1) Winkler, Gupta & Long (2002), (2) Slane et al. (2001), (3) Slane et al. (1999), (4) Ellison 
et al. (2001), (5) Allen, Petre & Gotthelf (2001), (6) This paper, (7) Green (2001), (8) Reynolds & Gilmore (1993) 
and (9) Combi, Romero & Bcnaglia (1999). 

''Galculated from expressions provided by Slane et al. (2001) and assuming a mean distance of 1.5 kpc, a volume 
filling factor for a sphere of one-quaxter and an explosion kinetic energy of 10^^ ergs. 



